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Abstract 

Neutrino oscillation experiments are known to be sensitive to Non-Standard Interactions (NSIs). 
We extend the NSI formalism to include one-loop effects. We discuss universal effects induced 
by corrections to the tree level W exchange, as well as non- universal effects that can arise from 
scalar charged current interactions. We show how the parameters that can be extracted from the 
experiments are obtained from various loop amplitudes, which include vertex corrections, wave 
function renormalizations, mass corrections as well as box diagrams. As an illustrative example, 
we discuss NSIs at one loop in the Minimal Supersymmetric Standard Model (MSSM) with generic 
lepton flavor violating sources in the soft sector. We argue that the size of one-loop NSIs can be 
large enough to be probed in future neutrino oscillation experiments. 
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I. INTRODUCTION 



There are many experimental results confirming that neutrinos have masses and oscillate 
between different flavors fT]. It is plausible that neutrinos acquire small masses from some 
high scale physics via the see-saw mechanism If such high energy dynamics are the 
only source of lepton flavor violation (LFV), neutrino oscillations could remain their only 
observable effects.^ It is possible, however, that extra sources of LFV are present at the weak 
scale. Such LFV sources could induce lepton number violating decays of charged leptons 
such as /i — )■ 67, /i — )■ eee and r — )■ /i7. Weak scale LFV interactions can also affect neutrino 
oscillation experiments. For instance, flavor violating Non-Standard Interactions (NSIs) can 
modify neutrino propagation in matter [SHTO]. 

NSIs could also affect the production and detection processes generating "wrong flavor" 
neutrinos [TTHT3] . Consider for instance an appearance experiment where a neutrino is 
produced in association with a muon, and a tau is detected at the detector. The standard 
interpretation of such a result is that it is due to — )■ u^- oscillations. This interpretation 
is correct as long as the produced neutrino is orthogonal to the detected one. However, if 
there are NSIs that violate the flavor symmetry at the source, the produced neutrino, which 
we denote by is not a simple flavor eigenstate. Analogously, in the presence of detector 
NSIs, the final state, which we denote as is not a simple flavor eigenstate. If these 
states are not orthogonal, that is, 

W)^e,^,^0, (1) 

tau appearance can occur without oscillation. 

In order to describe the effect of NSIs on neutrino oscillation experiments, let us discuss 
first a simple case with only two generations and a mixing angle 6 = tt/A. Consider also the 
case when the detector is at a distance L much smaller than one oscillation length, 

x^^«l. (2) 

In this case the oscillation amplitude simply reads Aosd^n ~^ ^t) ~ ix. When NSIs are 
present, they induce an extra contribution, ^nsi(^'m — >■ ^r) = ^^lt- Thus, to leading order in 
X, the total appearance probability is given by the squared of the sum of amplitudes 

P(Z/^ Ur) ^ \Aosc + ^NSlI^ ^X'^ + \e^,r? + Xm {Sf^r) ■ (3) 

The above simplified result has all the physics in it: the first term is the pure oscillation 
term, the second one is the x independent term that arises due to the NSIs, and the third 

^ A remarkable exception arises in the context of supersymmetric theories where such high scale dynamics 
could leave indelible footprints on the soft terms of the light sparticles via interactions not suppressed by 
inverse powers of the high scale [3 0] . 
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term is an interference term. Note that, when x ^ e^,-, the probabihty to detect a new 
physics effect is enhanced by the interference term, which is hnear in e^,-, compared to the 
typical quadratic dependence of lepton flavor violating decays of charged leptons [12] . This 
fact triggered renewed interest at present neutrino facilities [HHTQ] . The typical bounds 
on production and detection NSIs are about e^r ~ 10^^ [201 EI]- Higher sensitivities, of 
e^r ^ 10~^ are within the reach of future neutrino experiments [SSHSO]- 

In this paper we study NSIs at the loop level, extending the formalism of [HI [12]. We 
present a general framework that allows one to extract in a consistent way the physical 
parameters e^is (with a,P = 1,3) which arise at the loop level either from corrections to the 
tree level W exchange diagram or from more general corrections, in particular from scalar 
charged currents. We show how the physical parameter, e^^, can be obtained from the 
various loop amplitudes which include vertex corrections, wave function renormalizations, 
mass corrections as well as box diagrams. 

In the case of universal corrections to the W exchange amplitudes, NSIs emerge at one 
loop because after the Electroweak Symmetry Breaking (EWSB) the kinetic terms and the 
W couplings are generally not universal in the same basis. Rotating to the mass basis for 
the charged leptons, the misalignment between vertex and wave functions induce NSIs. We 
show that the associated Eajs are finite because the SU (2) l gauge symmetry protects them 
from possibly divergent contributions. 

This paper is organized as follows. In Sec. [Tl| we recall the standard formalism for NSIs. 
In Sec. |III[ we extend it to account for loop induced NSIs. In Sec. [TV] we discuss NSIs at 
one loop in the R-parity conserving Minimal Supersymmetric Standard Model (MSSM) with 
generic LFV sources in the soft sector, leaving the details of the calculation to the appendix. 
Finally, we present our conclusions in Sec. |V| 



II. NOTATIONS AND FORMALISM 

We start by setting our notations to follow Refs. [HI [I2]- We first consider only tree level 
processes and later discuss one loop effects. 

Neutrino mass eigenstates are denoted by i = 1,2,3, while |z/a) are the tree level 
weak interaction partners of the charged lepton mass eigenstates a~ , a = e, fi, r. These two 
bases are related by 

Wa) = ^UaiWi) . (4) 

i 

where U is the leptonic mixing matrix, the so-called PMNS matrix |31] . 

We consider experiments where neutrinos are produced at the source in conjunction 
with incoming negative or outgoing positive charged leptons. Then, the neutrinos travel to 
the detector where they are detected by producing negative charged leptons. We consider 
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new physics in the production and detection processes assuming that these NSIs have the 
same Dirac structure as the SM interactions and thus the amphtudes add coherently. We 
parameterize the new interactions at the source and at the detector by two sets of effective 
four-fermion couphngs, 

{GNp)a0 , {Gtsp)aP , (5) 

where a is the charged lepton index and /? is the flavor of the neutrino in the weak interaction 
basis. At tree level, the SU (2) l gauge symmetry implies that in the SM the four-fermion 
couplings are proportional to GpSa^. New interactions, however, allow for non-diagonal and 
non-universal couplings. 

Phenomenological constraints imply that the new interactions are suppressed with respect 
to the weak interactions. It is thus convenient to define small dimensionless quantities in 
the following way: 

^% = I p 2 P = s,d. (6) 

J\Gf + (G^p)aal^ + ^^^a l(^Np)a7l^ 

We denote by |z/^) the neutrino states that is produced at the source, and by the 
neutrino state that is detected 

At the leading order, we have 

el,-^^^, K) = M + ^:,b)- (8) 
The expression for the non-orthogonality parameter Sajs at the leading order is given by 

/ A + « = ,o^ 

where by C(e^) we refer to effects that are quadratic in or e^. Note that in the SM the 
non-orthogonality parameter vanishes, Ea^/s = 0. 

For simplicity, we consider now a two generation model where the production process is 
associated with a muon and the detection with a tau. We calculate the following expression 
for the transition probability 

P,r-MKm', (10) 

where z/*(t) is the time-evolved state that was purely at time t = 0. Using an explicit 
parameterization of the neutrino mixing matrix 

, cos 9 sin ^ \ 

^ = I (11) 
sin 9 cos 9 I 
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and keeping terms up to leading order in e we get 

P^r = sin^ X {sin^ 26 + 7^e(e^^ - e^^) sin 46} + sin 2x Xm(£^^) sin 26 , (12) 

where 



A 2 2 2 '-^ij-'-' /I o\ 

Am^^. = - m^., = x = Xi2 , (13) 



AE 

rrii are the neutrino masses, E is the neutrino energy and L is the distance between the 
source and the detector. 



A few remarks are in order regarding Eq. (12): 



We keep only terms up to 0{e). This is the reason why there is no effect at x = 0. 

The different x dependence of the various terms is important as it can be used to 
distinguish them experimentally. 

The interference term can be very important when x ^ e. 

The interference term depends on the imaginary part of the NSIs, that is, it requires 
CP violation. 

NSIs also affect the term proportional to sin^ x. Yet, within one experiment this change 
is absorbed into the definition of 6 and cannot be distinguished experimentally. 

In many cases the NSIs are closely related to lepton flavor violating charged lepton 
decays. However, they have a different dependence on e. Neutrino oscillation exper- 
iments are linear in e (if x ^ e) while decays like r — )■ fie~^e~ are quadratic. This 
makes neutrino oscillation experiments competitive in sensitivity. 



• With three generations the result is more complicated and can be found in |12j . 

Before concluding this section, we remark on the case of heavy neutrinos. For instance, 
consider the case of k heavy singlet neutrinos. The mixing matrix is 3 x (3 + fc) and the 
3x3 mixing matrix for the light neutrinos is not unitary anymore. In this case we have [32] 

h 

The point is that using neutrino oscillation experiments we can measure e, and claim de- 
tection of some new physics, but we can not disentangle the underlying mechanism which 
generates it. 

III. ONE LOOP NSI 

In this section, we extend the NSI formalism of Refs. [HI |T2] to include one-loop effects. 
In particular, we consider universal NSIs from correction to the tree level W interaction and 
non-universal effects due to box diagrams and scalar charged currents. 
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A. Correction to the W exchange amplitude 



At tree level, gauge invariance guarantees universality of the W interactions. This uni- 
versality is kept to all orders for an exact symmetry. For a broken symmetry, however, 
universality is lost beyond tree level. In the following, we show that one-loop effects make 
the W couplings and the kinetic terms of the fermions non-universal. In particular, we ex- 
plain why, in general, in the basis where the kinetic term are canonical, the W interactions 
are not flavor diagonal. 

In the following, we neglect neutrino masses since they give subleading effects to the 
NSIs, as we will discuss later. Similarly, we do not consider other possible non-universal 
flavor-diagonal NSIs. They can be there but they are assumed to be small and, to leading 
order, we can just add them to the effect we are considering here. 

On general grounds, one-particle irreducible one-loop effects include the self energy di- 
agrams for the charged leptons and the neutrinos, and the corrections to the W vertex as 
well. These one-loop diagrams modify the kinetic and mass terms for fermions and the W 
vertex by factors Z'^, Zj^^^, and 77^ defined as follows: 

- ^Pr-£,L7^(Zf)^V,^ + H.c. ^,J = 1,2,3, (15) 

where the symbol "o" refers to "bare" or unrenormalized mass matrices and Z"' are matrices 
defined as 

{Z%^5,j + irj% a^u,e,W. (16) 

Note that Hermiticity of the effective Lagrangian ensures that {r]^''^y — if-'^ while may not 
be Hermitian. While at tree level rf — 0, lepton flavor violating loop corrections introduce 
non-zero 77 's. In general, the off-diagonal elements of 77" are finite, while the diagonal terms 
diverge. Of course, the physics is finite and we discuss how the divergences cancel below. 
We define 

Z^L = ^\,vZ'l'^^,v ' -^i ^ ^IZ^R^i , (17) 

such that Zf^ ji are diagonal matrices of positive elements, and and are unitary matrices. 

We rescale the lepton fields to make their kinetic terms canonical 

_i _i _i 

^ {Zl) ' 4 ^ U {Zi) ' ^ (4) ' . (18) 

where (Z")~^/^ is shorthand for the diagonal matrix of element {Z"')'^^^'^ . The charged lepton 
mass terms become 

>C„.ass = -IjR (4) ' R\ {rn} + ^i) (li) + h-C (19) 
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The mass terms (19) can be diagonalized by two independent rotations 



where and Rm are unitary matrices. We obtain 



RL [Zr 



R\[ml + ^l]LAZ\^ 



(frii) 



(20) 



(21) 



where m£ is diagonal. After performing the rescahng (18) and the field rotation (20), the 



kinetic terms are canonical and the charged lepton mass matrix is diagonal, whereas the 
interaction terms are not diagonal. For later convenience, we rotate the neutrino fields as 
ul — )■ Lmi^L, in order to keep them as much aligned as possible with their charged partners. 
Note that this choice is allowed because when we study NSIs we can neglect neutrino masses. 
As a result, the VT-boson vertices become 



92 



92 



int 



where 



Ll Z 



L\ Z^ 



Zl 



(22) 



(23) 



Let us stress that eq. (23) is valid to all orders in perturbation theory. Finally, the relation 



to the physical parameter e^^ can be derived from eq. rt7| and is given by 



-.w 



(24) 



^{ZZ^\UZZ%p 

This formula has a simple interpretation. Zap/ \f[Z7^Yaa \vfi) is the normalized state pro- 
duced at the source, while its conjugate is the one detected. 

We proceed to find the leading order expressions for Z^^ and e^. That is, we will work 
to one loop level. In this case, we can identify the off-diagonal terms of Z with the NSIs at 
the source and the detector, and therefore we find 



— Cq,^ — c 

For the physical parameters we then get 

Zap + Zl^ 



d * 



'a/3 



(25) 



(26) 



The same result can be obtained directly from eq. (24). Note that ZZ'^ = dap when Z is 
unitary. When the deviation from unitarity is small, ZZ'^ = 1 -|- e^, we recover the previous 



result. See also eq. (33) below 



At one loop, the transformations for the lepton fields of eq. (|18|) read 
1 



i^L-> 1 



(27) 
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Similarly, to one loop accuracy, the transformations of eq. (20) read 

The unitarity of Lm and Rm implies that 

6LI^ = -5Lr. 
In this approximation, Z is given by 

Z = 1- 



SRl 



-SR„ 



Vl 



w 



(28) 



(29) 



(30) 



where rf , ?7^, and rf^ have to be evaluated at one loop accuracy. Eq. (|30|) shows that, at this 



level, the dependence of Z on drops completely out (its dependence would be reintro- 
duced at two loop). Then we learn that at leading order the non-orthogonality parameter 



between the source and detector neutrinos is given by (|26|) and reads 



JW 



s* , d 

tf. 



I w 
+ Vl 



where we used the fact that rj^ and ■i]'^ are Hermitian 



Vl 



Vli , 

ap 



(31) 



Eqs. (23), (26), and (31) are the main results of this section. A few remarks are in order 



when inspecting them: 



1. For any given model, eqs. (23 ), (26 ) show how to extract the physical NSI effects out of 



the calculations of the vertex corrections and the self energies encoded in the rotation 
mass matrices Lj^, and L^. However, eq. (31) demonstrates that, at one loop 
accuracy, all we need to do is to calculate t]'^, ?7^, and ?7^, since starts contributing 
only at two loop. 

We note that is finite. While each of the diagonal terms in r]"- may diverge, the 
combination is finite because of the SU (2) l gauge symmetry. This can be seen by 
the fact that the UV properties are insensitive to EWSB. Thus, the divergent part 
of 7]'^ is real, flavor universal and independent of a. Therefore, the divergent part 

of (,rt + 

vY - Vl~ vDaa, as well as its renormahzation scale dependence, cancel. 
In contrast, the flavor off-diagonal elements of r/]^ and r/^'^ are singularly finite and 
scale independent, as a result of the GIM mechanism. We stress also that when the 
electroweak symmetry is restored, vew ~^ 0, then the finite parts are flavor universal 
and e^^ — ?■ 0. Both results are illustrated in a concrete example below where the rj"' 
are explicitly calculated. 



3. Considering the CP-conjugated process, we obtain 



•^NSI 



-a/3 



(32) 



So, CP is violated when the rj"" have non trivial imaginary parts. 
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4. Eq. (31) admits an interpretation in terms of the scattering amplitudes Aiai3- The 
idea is to think about neutrinos as invisible intermediate states, and sum over all of 
them in the propagation process from the source to the detector. Considering a source 
neutrino associated with a charged lepton a and a detection that is done by a charged 
lepton (3, the scattering amplitude is given by 

Map oc (ZZt)^„ = + (riY + r/f t _ _ = + . (33) 

\ / pa 

where in the second step we expanded in rf/'^ ■ Thus, we learn that the non-universal 
part of the amplitude is just e^. 

B. Non-universal effects 

So far, we have considered only NSIs induced by the universal corrections to the W vertex 
and self energy diagrams. These effects are independent of the production or detection 
processes as long as these processes are mediated by W exchange. We now move to discuss 
other loop effects that are not universal, that is, that may be different for different production 
or detection processes. 

As mentioned before, the universal effects are SU (2) breaking effects and therefore 
suppressed like M'^/Mljp. In the effective field theory language, this would correspond 
to the effects induced, after the EWSB, by gauge invariant dimension-six operators like 
(Llt°-'^^Ll){H'^t°-D^H), where Ll (H) stands for the lepton (Higgs) doublet, are either 
the identity or the SU{2) generators and is the covariant derivative. 

It is then clear that contributions arising from dimension six four-fermion operators, 
also suppressed by M^/M^p, must be consistently included along with the W^-penguin 
contributions. Therefore, for any given New Physics (NP) model we can write the following 
expression for the physical Eais parameter 

^«/3 = ^5 + (4S"™"V(4J"-^ (34) 



where e^g has been already defined in eq. (31 ) and we have assumed both matter effects and 



higher order operators (with dim > 6) to be negligible. 

In this work we are interested in NSIs with the same Dirac structure as the SM in- 
teractions. The reason is that they maximally exploit the interference between SM and 
NP amplitudes. In particular, focusing on realistic production and detection processes like 
fi — )■ eUgiJa and P — > fiUa (with P = n,K), the relevant dimension six operators are the 
following 

^ + «J'^'"-') {^.1^PL^^) {nxPLVe) , (35) 



(5.a + {el^t"^-') {ui'PLd) {JUxPlVo) , (36) 



I (^d \dim-6 

V2 
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where (e^'^)"^*"*"^ stand for the loop-induced corrections. 

As we will discuss in the next section, in the context of supersymmetry might be 

induced either by means of gaugino/sfermion boxes or through the tree level exchange of 
charged Higgs with loop induced flavor changing couplings Hlv. 



Matter effects 



So far, we have discussed only NSIs induced by the charged currents. However, we would 
like to emphasize here that there exist also NSIs in matter via neutral currents, even for 
negligible SM matter effects. 

Indeed, starting from the SM neutral current interactions 
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eff 



2c 



w 



Z^-U^LY{Ziy u,L + }l.c. t, J = 1,2,3 



(37) 



where hereafter cw = cos^vy and sw = sm9w, after performing the rescaling (18) and the 



field rotation (20), the Z-boson vertex with neutrinos is given by 

Z = Ll{Z'r 



I LI Zl [Zl 



(38) 



where, again, we rotated the neutrino fields as Lm^i- At one loop accuracy, eq. (38) 

becomes 

Z = i+T^l-r^l. (39) 



We stress that also NSIs induced by the Z-penguin are SU(2) breaking effects and all 
the considerations made for l^-penguin induced NSIs apply also here. Therefore, besides 
Z-penguin effects, we have to include dimension six four-fermion operators. The relevant 
neutral interactions are 



2Q/S 



w 



m,f\dim—6 



(40) 



where / stands for the fermions in the matter. For normal matter, / could be electrons, 
protons and neutrons f = e,p,n, Qf is the electric charge of / and I^j^ is the third component 
of weak isospin of the left-chiral projection of /. 

Since we can identify e^^ = Z^^, the total physical parameter in the matter e^^-^ reads 



m,f\dim—6 



(41) 



When matter effects are included, the transition probability P^r of Eq. ( 12 ) becomes 

X sin 26' 



P ~ 



X 



sin 20 - 2L^A^ 7^e(e 



f 



+ x^ 7^e(e^^ - e^^) sin 40 + 2x Xm(£^^) sin 20 , 



(42) 
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where A-^ = y/2GFnf and we have assumed x <^ 1 and constant fermion densities nj. 

A close look to Eq. (42) shows that the interference term between the SM and non-SM 
matter effects (^p-^) depends only on the real part of e^^-^ . By contrast, for NSIs at the 
production or detection processes the interference term depends only on the imaginary 

part of the NSIs. 

The situation changes when considering the transition probabilities Pg/^ and Per which 
involve the electron neutrino. In these cases, there are interference terms, driven by charged 
current SM-effects, which are also sensitive to the real parts of e**'*^ 



D. Scalar charged current 

Many UV completions of the SM contain an extended Higgs sector, for example, the 
MSSM. On general grounds, the presence of at least two Higgs doublets leads to a mis- 
alignment between the fermion mass matrices and the Yukawa couplings. As a result, Higgs 
mediated FCNC processes are induced already at tree level resulting in large effects, unless 
a flavor protection mechanism is at work. Indeed, the Natural Flavor Conservation (NFC) 
hypothesis was introduced to deal with this flavor problem. However, even if NFC holds at 
the tree level, this hypothesis is spoiled by quantum corrections [33] . For instance, if NFC 
arises as a result of a continuous PQ symmetry, the breaking at the quantum level of such 
a symmetry (as it is required in order to prevent the appearance of a massless Goldstone 
boson) would reintroduce FCNC effects [53] . 

This is the case of supersymmetry where the holomorphy of the superpotential implies a 
type-II structure of the Higgs potential at the tree level. Yet, the presence of a non- vanishing 
//-term (such that fiHuHd) induces, after SUSY breaking, non-holomorphic Yukawa cou- 
plings for fermions (such as QiduHl) [34j and therefore Higgs-mediated flavor violation is 
unavoidable [35] . 

Bearing in mind the above considerations, in the following we perform a model indepen- 
dent analysis of NSIs arising from loop-induced scalar charged currents. 

The charged Higgs couplings with leptons are described by the following effective 
Lagrangian 

Aff = y-jL {Ve + V^^y ^^RH^ + H.c. (43) 

where 

y°,= mnan/3 (44) 

and we recall that "o" refers to unrenormalized quantities. 



The transformations of eqs. (18), (20) lead to the following effective Lagrangian valid to 
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all orders in perturbation theory 



-\ 31' 



H+ + H.c. 



(45) 



In order to find the leading expansion for the above effective couplings, we proceed as follows. 



We first rescale the lepton fields at one loop level, see eq. (18), so that the charged lepton 
mass terms become 



where 



^mass = -^jR {m} + (Jm^)^* In + h.c. , 



(46) 



(47) 



Then, we make use of the one loop expansions for the matrices and Rm of eq. (28) 
leading to 



1 J* 



+ 6me — 5Rmm°f + m\5Lm 



m 



I 1 



(48) 



where is diagonal and we have consistently retained only the leading one-loop terms. At 



this level, the unitarity condition of eq. (29) ensures that Re(5L:^) = Re(5-R:^) = and the 
corrected mass eigenvalues are given by 



mi. = m°i. + Re {,5mi,).. , 

while the condition of reality for the masses implies 

Im {Smi)jj 



Im 5R" - Im 6L" 



ma 



(49) 



(50) 



Finally, SL^ and SRm are determined by 



5U1 



mt. — mi 



(51) 



SRli 



miASmly + {dm^y^mi^ 



m. 



m-n 



(52) 



where m^. = {m£)ii, that is, it is the ith eigenvalue of me. We are ready now to expand 
eq. (45) up to one loop. By making use of the eqs. (47), (49), and (50), we obtain the 



following flavor conserving couplings 



1 1 Tl^t 

^ ' ^ 2 2 mt. 



+ r] 



H 



,H+ + H.c. 



(53) 
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where 



tan (3. 



(54) 



w 



Again, while each term in eq. (53) is in general divergent and renormalization scale depen- 
dent, their sum is finite and scale independent. 

For the fiavor violating charged Higgs couplings, we find the one loop expression 

1 ji 



H+ 

eflf 



n 



■ V 1 I 



H 



iy+ + H.C. 



(55) 



and each term in eq. (55) is finite and scale independent thanks to the GIM mechanism. 
Notice that, in contrast to the case of NSIs at the W-hoson vertex, we are now sensitive to 
the rotation mass matrices Lm and Rm already at the one loop level. As we will discuss in the 
next section, within the MSSM the rotations 5Lm and SRm actually provide the dominant 
effects to the flavor changing couplings. 

Let us consider now the case of j = 3, that is relevant for a tau neutrino production. In 
such a case, the one loop expansions for 6Lm and 6Rm of eqs. (51) and (52) take the form 

3i 



3i 



{5Rm, 



m^ 



1 e 
+ -/r 

+ 2^i 



rriT 



+ — % 



3i 



Finally, inserting the above expressions for bLm and bRm in eq. (55), we flnd 

= I7,i Zl + H.c. , 

where 



ry'ii 



1 1 rf^ 



3i 



(56) 
(57) 

(58) 
(59) 



Let us remark that NSI effects driven by charged scalar currents are expected to be par- 
ticularly relevant for the neutrino production via charged meson decays. In fact, whenever 
the relevant Yukawa couplings are proportional to the fermion masses, only processes like 
P iu (with P = K,7t), which are helicity suppressed in the SM, might receive large con- 
tributions. Other production processes like fi — )■ ez/z/ or detection cross-sections are expected 
not to be signiflcantly affected by such charged scalar currents. As a result, we now have 



and therefore e^r ~ e^*- This is in contrast to the case with dominant NSIs at 



-a/3- 



the W-boson vertex where, as we already discussed, it turns out that e^^ 

We can proceed now to establish the relation between e^,- and in the case where the 
neutrino source is given by the process tt ^ fiu. This decay is mediated by tree level 
and exchanges. The relevant effective Lagrangian describing this process is 



4Gf 

71 



VUui^PLd){jirPLi^^) + V, 



ud 



ydZ_ 



32* 



m 



uPRd){flPLl^r) 



(60) 
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where Pr^l = (1 i 75)/2 and yd is the down quark Yukawa couphng. Since the tt meson 
is a pseudoscalar, its decay amphtude can be induced only by the axial-vector part of the 
coupling and the pseudoscalar part of the coupling. Then, once we implement the 
PCAC relations 

///" 

(0|M7/i75C?|7r") = i/^p^^ , (0|M75(i|7r") = -if^ -f — , (61) 



we get the amplitudes 



G 

MZ^^,^ = -^VudUm^ Ji(l - 75)i^/. , (62) 



We observe that the SM amplitude depends on the lepton mass because of the helicity 
suppression, in contrast to the charged Higgs amplitude which does not suffer from this 
suppression. Yet, many NP models predict the Yukawa couplings to be proportional to 
the fermion masses. Effectively then, in such models both the PF-boson and charged Higgs 
amplitudes have the same lepton mass dependence. 

Finally, recalling that the produced neutrino state is = + e^^|i/r), we identify 

In the case of K ^ fiu, the relevant parameter e^^ can be simply obtained from through 
the replacement {mT^,y(i,md) {ttt-KjUst^s) and we get 

^It ml ms + rriu ya 1 

~1F — — I ~ ^ ■ (65) 

mj^ md + rriu ys 20 

It has yet to be seen which process tt ^ iii/ or K ^ /ii/ may represent the best probe of this 
scenario when combining the more intense neutrino flux obtainable from tt — > //i/ with the 
higher NP sensitivity ol K ^ jiv. 



IV. ONE LOOP NSIS AND SUPERSYMMETRY 



In this section, we apply the model-independent formalism developed in the previous 
section to the case of the R-parity conserving MSSM with new sources of LEV in the soft 
sector. We will analyse first loop-induced NSIs from the V — A charged current, passing then 
to NSIs from the charged scalar current induced by the heavy Higgs sector of the MSSM. 
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A. V — A charged current 



In the MSSM NSIs may be induced by the V — A charged current through W^-penguin 
as well as box contributions. The former effect arises only after the EWSB and therefore is 
suppressed by M^/Mj^p. In particular, within the MSSM, there are three possible sources 
of SU{2) breaking: 

i) the D-terms, 

ii) the left-right mixing terms, and 

iii) the neutralino/chargino mixing terms. 

The latter comes from dimension six four-fermion operators and is also suppressed by 

The full analytical calculation relevant for NSIs in the MSSM, upon which our numerical 
analysis is based on, is reported in the appendix. In the following, instead, we prefer to 
discuss general properties within an illustrative toy model which is a particular limit of the 
MSSM that greatly simplifies the calculation but still retains the most relevant features. We 
use standard MSSM notation (for a review see for example, Ref. [36]) and we consider only 
the lepton sector. 

In our toy model, we decouple the Higgsinos, H, by taking a large yU-parameter and we 
decouple the Winos, W, by taking their soft mass term, M2, to be very large. We also 
assume left-right mixing between the sleptons to be negligible. Finally, we take the EW 
vacuum expectation value, vew^ small compared to the soft mass term of the Bino, Mi, and 
the soft term for the left-handed sleptons. Ml. That is, we consider a model where 

A m^/itan/3 ^ ,^ ,^ 

TT ' ^,2 < 1 ' VEW-^Mi^ ML'^^^^ M2. (66) 

where A stands for the trilinear soft terms. In this model there is only one neutralino, the 
Bino. In practice, the model looks supersymmetric only with respect to f/(l)y. The EWSB 
can be treated as a perturbation. Note that, in our toy model, only the SU{2) breaking 
source i) is at work. Later on, we will also discuss the impact on NSIs of sources ii) and iii), 
which are expected in a more realistic model. 

The relevant interactions are the lepton-Bino-slepton vertex and the slepton-slepton-PF 
vertex, that are given by 

Cnsi = -V2ig'qYx'' {^lU^u'^i + ilUfl?) - V^ig {w^d^pd*^) + h.c. (67) 

where qy = —1/2 is the hypercharge of the left-handed leptons. U^^ (f^|*) is the unitary 
matrix that diagonalizes the sneutrino (charged slepton) mass matrix. In our model, the soft 



15 



FIG. 1: 1-loop contributions to the lepton self-energies (left) and to the vertex (right). The virtual 
particles running in the loop are sleptons and a neutralino. 



terms are SU{2)l x U{1)y symmetric and, since there are no left-right mixings, Ua = Ui- 
Thus, from this point on we use Ui for both terms. 



According to Eq. (31), all we need to calculate are the loops in Fig. [ijand extract t]^, rf 
and . In our calculations we use a naive UV cutoff, that is, we perform the /c^ integral 
up to A^. We further introduce an unphysical mass parameter [i to make the arguments of 
all logarithms dimensionless. Effectively, this is equivalent to introducing a renormalization 
scale. As a check on our calculation we see that the final physical results, that is, e, is 
independent of these two unphysical parameters. 

We first calculate the wave function corrections. The only diagrams contributing to ri" 
at one loop are through the exchange of the Bino and sneutrinos. The result is 



ri 



V 

k ) 



(68) 



where 



167r2 
-log 



A2 / 1 



2 



A2 



+ 



4 



(m 



m 



2 \2 



log 



I^k 



ml + m^„ 
2(m? — m^o) 



(69) 
(70) 



The calculation of the wave function correction to the left-handed charged leptons rj^j^ 
proceeds in a completely analogous way, the only difference being that now we have sleptons, 
instead of sneutrinos, running in the loop. Therefore, r]j^ reads 



(71) 



where is simply obtained from by the replacement m?^ — )■ m| . 

Next, we calculate the one loop correction to the W vertex, rf^ . The result is 



w 



/2 2 

9 <1y 



kj*jjki j-W 



k 5 



(72) 
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where 



(73) 



1 



i2 



m 



<(mJ.-,4)lo; 



+'^iS'^lk - "^?») log ( ^ ) + <o(m| - m^J log 




(74) 



In order to obtain the final result we use Eq. rt31| with Eqs. (68), (71), and (72). We get 



k 



where we neglected (9(A^^) effects. 

We are now in a position to check the finiteness of the physical amplitude. Note that 77'", 
?7^, and rj^ contain a log-divergence. SU{2)l gauge symmetry constraints the coefficients of 
these two divergences to be equal to each other. We can see that this is indeed the case by 



direct inspection, e in Eq. (75) depends only on the functions that are independent of A. 



We can also check that the result is independent of fi, as it should be. For this note that the 



sum, 2F^- 



, is independent of /i. While the above results are automatically achieved 



by each off-diagonal term contributing to e, as a result of the GIM-mechanism, their validity 
for the diagonal components represents a check of the correctness of the calculation. 

Another important check is to make sure that in the limit of no EWSB, there is no effect 
induced by the kinetic term because SU (2) ^ gauge symmetry makes it aligned with the 
ly-interaction. When sending vew ~^ 0, the charged sleptons become degenerate with the 
sneutrinos, m? 
to the identity 



m^^. In this limit, using (|75|) we learn that the relevant sum is proportional 



/ W , Wi V i \ 



OC UjU^ 



6, 



al3 ■ 



(76) 



vew=0 

We learn that no ffavor changing amplitude is induced thanks to the unitarity of f/^. 

The fact that the effect vanishes for vew = can be used to get an approximate formula. 
We can define the presumably small parameter 



that vanishes for vew = 0, and expand in to the leading order 



W I u t 

7] +7] —7] — ?7 



167r2 



(77) 



(78) 
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where Gk is a function of SUSY masses which does not vanish in the hmit of Ofe — > 0. The 
dominant sphtting between left handed sneutrinos and sleptons originates from the D-terms 
and is flavor universal 

(m?^ -mjj^ ml cos^ Ow cos(2/3) . (79) 
As a result, in this toy model, e^r can be estimated as 



where we have defined the mass-insertion parameters 5^, = {M^ij/ (MDaiMDjj, as usual. 



Such values are well below the expected experimental resolutions even for ~ 1. 



We discuss now the main differences arising in the cases where the SU (2) breaking sources 
ii) and iii) are switched on. In the case ii), the leading effect for reads 

ay ml \Ar - ^itan 13 f L /qin 
4:71 m^^omj ^ 

where we picked up a double left-right mixing term for the third slepton generation. In 
principle, e^r could reach values even slightly above 10~^ for m^// tan /3/m| ~ 1; in practice 
the constraint from r — >■ //7 implies that Sfj^r < 10~^. 
Finally, in the case iii) we get 

where are the mixing angles of the chargino mass matrix which read 

where tan/3 = v^/vd- We have explicitly checked that also in this case e^^ < 10~^ after 
imposing the constraint from r — >■ /i7. 

We discuss now the box induced NSIs. These effects survive even in the hmit where 
all the SU{2) breaking sources discussed above are negligible. In particular, it turns out 
that the largest effects arise for light sleptons/ Winos and heavy Higgsino/Bino. The latter 
condition is necessary to suppress BR(t — >■ /i7). As a result, it turns out that 

As we will show in the numerical analysis, the box contribution provides the dominant 
effect to Efj^r that can reach experimentally interesting levels while being still compatible 
with the current bound on BR(r fij). In fact, in the most favorable situation where 
M2 ~ mi <^ II Ml, one can find that 



as we will confirm numerically. 
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B. Scalar charged current 



In the MSSM, Higgs mediated LFV effects are generated at the loop level, e.g. see 
Fig. |2} In fact, given a source of non-holomorphic couplings, and LFV among the sleptons, 
Higgs-mediated LFV is unavoidable 




FIG. 2: A contribution to the effective UrinH^ couphng. 
Starting from the model-independent parameterization for the flavor violating couplings 



of a charged Higgs with leptons, eq. (55), we specialize now to a SUSY scenario assuming 



R-parity conservation to avoid tree level flavor changing effects. 

Analysing the full expressions of such couplings (reported in the Appendix), we find 
that the field rotations 6Lm and 6Rm. induce the dominant contributions for the effective 



Lagrangian of eq. (55). In fact, one can show that they are parametrically enhanced by a 



tan/3 factor 6Lm,SRm ~ [«2/4vr] x tan/3 compared to ?7^'^ ~ a2/4vr and t]^ ye x [0:2/47?]. 

Since the effects we are going to discuss can be relevant only if tan/3 ^ 1, it turns out 
that 

C^s - {y\ - ^L^yl + yl^Rn.f i^RH+ + H.c. . (86) 

Retaining only the dominant tan /3 enhanced contributions in the corrections to the lepton 
mass matrix, one has that 6mi ~ r/,^ and therefore 



m 



(87) 



where e,eL,R are loop factors of order a2/4vr which depend on SUSY mass ratios and tf^ 



tan /3. Therefore, the rotation matrices can be determined explicitly from eqs. (51 ) and (52) 
and they read 



6L 



3i 



6R 



3i 



(88) 



:i + etfs) ^ ' (1 + et/s) - ■ (1 + et^) 

Finally, in the basis where ul — > Lmi^i, the effective Lagrangian for the if^ couplings with 
leptons reads 

92 



(89) 
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where we have defined ^''l(r) = ^l{R)^l{r)/(^ + 

An inspection of the above effective Lagrangian reveals that: 1) since the Yukawa operator 
is of dimension four, the quantities A-^*^ depend only on ratios of soft SUSY masses, hence 
avoiding SUSY decoupling. Yet, the NP effects induced in physical observables will decouple 
with the charged Higgs mass; 2) the loop induced flavor violating couplings are enhanced 
by an extra tj^ factor compared to the tree level flavor conserving couplings. Therefore the 
loop suppression can be partially compensated if A-^*^^^t^ ~ 1; and 3) the flavor violating 
couplings H'^ii'r (with i = e,fi) exhibit a Yukawa enhancement factor mr/me compared to 
flavor conserving couplings H^iui when they are induced by A-'* 



R- 



Applying the above results to the model-independent parameterization of eq. (64), we 
flnd 

where is a non-holomorphic threshold correction stemming from the quark sector typically 
of order ~ 10^^. 



As seen in eq. (65), it turns out that e^^/e^ ~ 1/20. We notice that and show an 
enhanced sensitivity to sources of flavor violation in the right-handed slepton sector thanks 
to the Yukawa enhancement factor 771^-/771^. 



In order to quantify the allowed size for e^^^ , we have to impose the constraints arising 
from the charged lepton LFV decays. The most sensitive probe of Higgs mediated effects is 
generally t ^ ij7] [38] and the corresponding branching ratio reads 



^32 


2 + 


Af 







10-^ ' V^' H , . i-Tm , . I. > (91) 



Br{T ^ IJU^Ur) \ 77l\ J \l + tqtp\'^\l + tt 

where 771^ is the pseudoscalar mass such that 7n\ = — at tree level. Imposing the 
experimental constraints from i?r(r — t- fiTj) < 10^'', it turns out that 

< 10-2 , e;^<5x 10-^ , (92) 

where the above bounds arise for IA32I ^ |A^^|. 

Finally, let us mention that Higgs mediated LFV interactions also induce lepton univer- 
sality breaking effects in P ^ iu {i = e, /i) [39] . However, these effects can only constrain 
|A|j I which is unrelated, in general, with the relevant LFV term for NSIs, that is |A|j^|. 



C. Numerical analysis 

In this section, we provide the predictions for the NSI parameter in the framework of 
the R-parity conserving MSSM with generic LFV soft breaking terms. The allowed values for 
X77i{e^^r) are obtained after imposing the following constraints: i) the data on flavor physics 
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observables; ii) the mass bounds from direct SUSY searches; iii) the requirement of a neutral 
hghtest SUSY particle; iv) the requirement of correct electroweak symmetry breaking and 
vacuum stability; and v) the constraints from electroweak precision observables. 

Concerning NSI effects driven by the charged Higgs exchange, the most stringent bounds 
come from the data on LFV and i?-physics observables. In particular, the processes B 
Xgj, B ^ TV and B — t- Dtv are known to be the most powerful probes of new charged 
scalar currents. In principle, also the process Bs^d ^~ shows an enhanced sensitivity to 
extended Higgs sectors. However, since the loop-induced flavor changing coupling Hbs{d) 
(with H = H^, A^) depends on the details of the soft sector, to be conservative, we do not 
impose here the (model-dependent) constraint from BR^Bg d — )• /i^/i"). 

The bounds from BR(i? — )■ Xg'y) have been obtained employing the SM prediction at the 
NNLO of Ref. gO], BR{B X^r, > 1-6 GeV)™ = (3.15 ± 0.23) x 10-^ combined with 
the experimental average |1IH13] BR(5 ^ X,r, > 1-6 GeV))"^P = (3.55 ± 0.24) x 10"^ 
As for the SUSY contributions, we use the calculation of Ref. [Sj assuming decoupled 
gluinos and squarks. For B — rz/, we use the current world average BR(i? — )■ rz/)exp = 
(1.73 ±0.35) X 10"^ t45j, the SM prediction (1.10 ±0.29) x 10"^ [^6] (see also |17]) and the 
NP contributions of Ref. [IH] . 

Finally, the NP sensitivity of B Dtu can be better exploited normalizing it to BR(i? 
Dtu)/BR{B Diu) where i = e, fx We use the world average (49 ± 10)% [SI] and 

the theoretical prediction of Ref. . 

In our numerical analysis we impose all the above constraints at the 2a C.L.. 

In Fig. ^ on the left, we show the values attained by |Ime^|, see eq. (90), in the 
tan/3 — Mh± plane setting the LFV parameter |A^^| = 10^'^ (varying A^^, |Ime^| would 
rescale according to |Ajj^|/10^'^). The red, green, glue and yellow regions are excluded by 
the current bounds on i? — )■ Xg'-f, B — t- tu, and B — t- Dru and r — t- fir], respectively. 

As shown by fig. (pi), |Im£:^| can vary in the range (10~^, 10~^) for tan /3 < 60 and Mj^i < 



500 GeV. The corresponding values for jlmej^^l can be obtained by |Im£:^^|/|Im£:^| ^ 1/20. 

We now discuss the NSIs as induced by the V — A charged current via the one loop 
exchange of gauginos/sleptons. As discussed in the above section, the dominant effect to 
liae^r arises from the box contributions. In fig. ^ on the right, we show the correlation 
between BR(r — t- fi'-f) vs. jlme^^l in the case where the neutrino source is provided by the 
muon decay fi — )■ euri^e- We have assumed heavy squarks implying negligible NSIs at the 
detector level. In this limit also NSIs for the production process P — ?■ are suppressed. 
Moreover, as the largest effects for Ime^r are obtained for light sleptons/Winos and heavy 
Higgsino/Bino (to keep under control BR(r — /i7)), we employ the following scan over the 
SUSY input parameters: M2, mi < 1 TeV, fi, Mi > 500 GeV and 3 < tan /3 < 10. 

As shown by fig. (|3]), llme^^-l can reach experimentally interesting values |Im£:^T-| < 3 — 
4 X 10~^ and this would unambiguously imply a lower bound for BR(r — )• fi'y) quite close to 
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FIG. 3: Left: NSIs in the process K ^ induced by Higgs mediated effects. Right: NSIs in the 
process fi — )• ez^z/ induced by VF-penguin and gaugino/slepton boxes. See the text for details. 

the current bound. 

V. DISCUSSION AND CONCLUSION 

The idea that neutrino oscillations can probe NSIs is very attractive. In theory such 
experiments are sensitive to any form of new physics that makes the produced and detected 
neutrinos non-orthogonal. Such non-orthogonality, parameterized by Sa/s, niay come from 
new tree level interactions, new heavy neutrinos, or one loop effects that modify the couplings 
of the W boson to the leptons. 

In this work we presented a general framework that allows one to extract in a consistent 
way the physical e arising at the loop level either from the V — A or scalar charged currents. 
We show how e can be obtained from the various loop amplitudes which include vertex 
corrections, wave function renormalizations, mass corrections as well as box diagrams. 

As an illustrative example, we discussed NSIs in the R-parity conserving MSSM with new 
LFV sources in the soft sector. 

We argued that, in general, the size of one-loop NSIs is quite small, e ~ O{10~^). To be 
observed, such small numbers require very precise measurements of the neutrino appearance 
probability as a function of L/E. We hope that such measurements will be possible in the 
next generation of neutrino oscillation experiments. 
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VI. APPENDIX 

In the following, we provide the full analytical expressions for the self-energies and vertex 
corrections relevant for NSIs in the R-parity conserving MSSM. For the Feynman rules, we 
closely follow the notation of Ref. [53] . 

The lepton self-energies read 



{VvrY = RifMRiLNBi{mL,.m^^) + R^R^ B,{m,,,mc^ (95) 
(47r)2 [riUy = -Li'r^RitNBoimL,,m,,J - L'^cRI'c Bo{m,„mc,) . (96) 



IJ 



LiLNLiLNB^{mL„mr,J + Li^lL'Jl^Biim^.^mcJ 



(94) 



The vertex corrections relevant for Wiu are 



+ 



+ 



+ 




(97) 



The vertex corrections relevant for H^v are 




- R^NCH^Cr^^N, Co {nil^ , ml^ , ml 
+ Lif^Rilc [L%HC2{ml,m%,mlJ 

- R^NCH^Ni mc^ Co{ml^, m% , m^^ 



)] 



)] 



(98) 
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The gaugino/slepton box diagrams relevant for the process /i — > ei^rfe read 

The box diagrams generated by the gaugino/slepton(squark) exchange contributing to 
production process P — > ixv^ {P — tt, K) and to the detection process read 

+ L/^-L^J^LS-W*^d&c^2(m,,, m^,, mc„ m^r,) 
The expressions for the loop functions appearing in the above amplitudes read 



So(mi,m2) = - + 1 
e 



rrin 



2 1 ^1 2 1 
iog^ — 1712 



Bi{mi,m2) 



+ 1 



1 "^2 
log^ + 



ml 



2 2 
mf — 1712 



mi 1 mf + mi 

log^ + 1 

mf 2 mf — m2 



Co{mi,m2,m3) = 



m? — ms 



mo 



mt 



ms 



log 



ms 



m? 



mi 



mt 



ms 



log 



m? 



mr 



mg log ml /ml 



^ , , 1 ^ , m? mologml/m? 

C2(mi,m2,m3) = - + 1 + log — + 2n7-^ 2V + ^ 2 2n^ 2 2n ' 

£ //^ [mf — m2){m^ — m2) [mf — m^)[m2 — m^) 



Do(mi,m2,m3,m4) = 



mf log mf 



(m^— mf)(m|— mf)(m|— mf) 
+ {lo2}+{lo3} + {lo4}, 



D2{mi,m2,m3,m4) = 



1 



mf log mf 



4 (m| — mf ) (m| — mf ) (m^ — mf ) 
+ {lo2}+{lo3} + {lo4}. 
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